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In [lj], Rashedi, Nezamabadi-pour and Saryazdi propose an interesting 
optimization algorithm called 'GSA: A Gravitational Search Algorithm'. The 
paper has received numerous citations since its publication in 2009. 

GSA simulates a set of agents in an iV-dimensional search space, which 
behave as point masses. The position of agent i, denoted by Xi, represents a 
candidate solution for the optimization problem at hand. The mass of agent 
i, denoted by m^, is correlated with the quality of the candidate solution 
it represents: the higher the quality, the higher the mass. The agents are 
said to "interact with each other based on the Newtonian gravity and the 
laws of motion" [lj. If the law of gravity were to be strictly followed, the 
magnitude of the force exerted on agent i by agent j (and vice versa) should 
be "proportional to the mass of each and [vary] inversely as the square of the 
distance between them" [2J. Formally, 



where G is a proportionality constant analogous to the universal gravitational 
constant and Rij = 1 1 X j oc 2 1 1 IS the Euclidean distance between agents i and 
j- 

The authors claim to have implemented a law of gravity where the magni- 
tude of the force between two agents is inversely proportional to the distance 
R between them. They use U R instead of R 2 , because according to [their] 
experiment results, R provides better results than R 2 in all experimental 
cases" [lj]. 
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In GSA, according to Equation 7 in [lj, the force exerted on agent i by 
agent j along dimension d G {1,2,..., N} of the search space is calculated 
as follows: 

where e is a small positive constant used to avoid division by zero. 

However, in order to obtain the force vector /y, its magnitude, 
should be multiplied by the unit vector pointing from agent i to agent j. 
Equation [2] is not correct in that (Xj — xfj is not normalized. The denomina- 
tor of Equation [2] should read i?? in order to implement an inverse-square law, 
or Rfj in order to implement an inverse-linear law (as the authors intended). 

Using Rij effectively makes the magnitude of the force independent of the 
distance between agents i and j. Formally, assuming e = 0, 

\\fij\ \ = Gniimj. (3) 

We do not in any way intend to undermine GSA's capabilities as an opti- 
mization algorithm; rather, our observation is that it is not strictly motivated 
by how gravity works (compare Equations [1] and [3]) . Mass and distance are 
both integral parts of the law of gravity; in contrast, the force model used in 
GSA is independent of the distance between the agents and is based solely 
on their masses (qualities of candidate solutions). 
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